Multi-objective pathfinder algorithm for multi-objective optimal power flow problem with random renewable energy sources: wind, photovoltaic and tidal

In this paper, the multi-objective optimal power flow (MOOPF) problem optimization objectives focus on four optimization objectives: generation cost, emission, real power loss and voltage deviation (VD). Three renewable energy sources with successful industrial applications, including wind energy, solar energy, and tidal energy are introduced. Renewable energy supply is uncertain, so Weibull distribution probability, lognormal probability and Gumbel probability are used to calculate the instability and intermittency of wind energy, solar energy and tidal energy, respectively. The inclusion of four energy supplies on the IEEE-30 test system and the consideration of renewable energy reserves and penalty cost calculation improve the realism of the model. In order to obtain the control parameters that minimize the four optimization objectives, a named multi-objective pathfinder algorithm (MOPFA) based on elite dominance and crowding distance was proposed to solve this multi-objective optimization problem. Simulation results show the feasibility of the model, and MOPFA can get more evenly distributed Pareto front and provide more diverse solutions. A compromise solution was selected by the fuzzy decision system. Comparison with the recently published literature also shows that the proposed model can effectively reduce emissions and other indicators. In addition, the statistical test results show that MOPFA's multi-objective optimization performance ranks first. In solving this complex optimization problem, results show the MOPFA is superior to other multi-objective algorithms in optimization accuracy and speed.

The stability of the electric power system, which serves people's production and lives, is an important problem. It is well worth studying how to optimize the parameters in electric control system. Optimal Power Flow (OPF) is an important tool for optimizing the power system, which is important for the reliable operation and cost reduction of the power system 1 . The OPF problem is characterized by nonlinearity and multiple constraints, these constraints include generator capability, line capacity, bus voltage and power flow balance 2 . The OPF problem belongs to the NP-hard problem, which adds to the difficulty of searching for the solution 3 . The goal of optimization is to find the optimal solution that minimizes the objective in the solution set of many control variables to be selected. Researchers usually consider the optimization of a single objective, such as the total fuel cost of generators, power loss, emissions and other objectives. However, reducing power loss should also reduce generator costs or other cost targets, therefore, multiple optimization objectives should be considered in the OPF problem, but unfortunately these goals are in conflict. The OPF problem of using thermal generators for power output has long been considered. On the other hands, with the widespread application of renewable energy, the

Related work
Traditional algorithms for solving the OPF problem. Since the OPF problem was raised 7 , many researchers have studied the OPF problem, and traditional methods for solving it include the Newton method 8 , quadratic programming 9 , linear programming 10 , and interior point method 11 . These methods mentioned in the previous question were used to solve the Optimal Power Flow problem in the early days, but these methods can only be solved for linear objective functions, which in turn can lead to an increase in error. In considering nonconvex fuel costs with threshold effects, these methods will not be solved efficiently.
Metaheuristic algorithms methods for solving the OPF problem. Researchers now have new solution ideas thanks to the emergence of metaheuristic algorithms, which have the benefit of not requiring them to concentrate on the objective function and constraints. Metaheuristic algorithms are widely used in solving OPF problems. In 12 a new genetic algorithm for coding systems was applied to the OPF problem, the objective was to minimize fuel costs. The particle swarm optimization (PSO) 13 algorithm was used to solve the OPF problem and test it on the IEEE 30-bus. In 14 , three new particle swarm optimization algorithms were used to find the optimal steady-state performance of power systems. Reference 15 provided a new initialization method for the problem of genetic algorithms that may be ineffective if starting values of voltage angles are selected quite randomly. Mahadevan et al. 16 applied a method named comprehensive learning particle swarm optimization (CLPSO) to the OPF problem with active power losses as objective functions. A new hybrid algorithm is proposed in 17 optimal reactive power dispatch problem with discrete and continuous control variables. In 18 , a Gaussian bare-bones water cycle algorithm (NGBWCA) was proposed to minimize resistive losses and voltage deviations. In 19 , authors solved the Optimal Power Flow problem using a modified Sine-Cosine algorithm, this algorithm obtains a lower value of fuel cost and power losses. In 20 , Wei and Zhou et al. employed an improved slime mold algorithm (ISMA) to find optimal control parameters in power systems, and its effectiveness and robustness were also demonstrated. Recent literature 1 proposed a high performance crisscross search based grey wolf optimizer (CS-GWO) to solve the OPF problem, fuel cost with valve-point effects and basic fuel cost are considered separately. Although the above study achieved the expected economic benefits, many researchers are not satisfied with achieving one optimization objective, and many researchers have also investigated the OPF problem with multiple optimization objectives. In 21 , multi-objective adaptive immune algorithm (MOAIA) was proposed for optimal reactive power flow incorporating static voltage stability. Reference 22 provided a new variant of the differential evolutionary algorithm, ensures high convergence speed and diversity of Pareto solutions, www.nature.com/scientificreports/ and extracts the best compromise based on fuzzy set theory. Pulluri et al. 23 proposed ESDE-MC methods to solve multi-objective OPF problem, non-dominated sorting and crowding distance was used in this method, the objectives to be optimized for the conflict include fuel costs, emissions, L-index and power losses. An improved NSGA-III (I-NSGA-III) was developed in 24 to solve multi-objective OPF problems, the optimization goal is to simultaneously minimize total fuel cost, total emissions, voltage magnitude deviation and power loss. Although many classical meta-heuristic algorithms have been applied to solve this problem, many excellent algorithms have been proposed in recent years, such as the cheetah optimizer 25 , etc. Therefore, the use of novel algorithms to effectively improve the accuracy of problem solving is still worth investigating. In recently published literature 26,27 , the improved heap optimization algorithm and multi-objective Manta Ray Foraging Optimization were developed separately based on the Pareto concept, these two algorithms aim to simultaneously optimize four metrics: fuel cost, emissions, power loss and voltage deviation.
A summary of past research shows that many researchers have conducted exhaustive studies on singleobjective and multi-objective problems for the OPF problem, and have obtained good results, but power systems are becoming more complex with the grid integration of renewable energy sources. Renewable energy sources have an element of uncertainty, yet the non-polluting nature of these sources compels us to use them, therefore, it is increasingly important to study the OPF problem with renewable energy. In 28 , the authors used the Weibull distribution to simulate the variability of wind, and then proposed a modified bacterial foraging algorithm to solve the OPF problem. On the basis of integrating wind power into the power grid 29 , also adding solar energy to the power grid, Lognormal probability distribution functions were used for forecasting solar photovoltaic power output. Reference 30 considers the OPF problem with or without wind power and tested it on IEEE 30-bus, IEEE 57-bus and IEEE 118-bus respectively. In 31 , a novel hybrid modified imperialist competitive algorithm and sequential quadratic programming were proposed to solve the OPF problem, which studied the uncertainties of solar power and wind energy. Li and Gong 32 proposed an enhanced adaptive different evolution and then applied it to the OPF problem on a modified IEEE 30-bus system, in which wind power and photovoltaic power are also being considered as energy supplies. A new version of the JAYA algorithm was proposed 33 to solve the problem of OPF incorporating renewable energy sources, using a linear weighting method to integrating a multi objective OPF problem with four different objective functions into a single objective optimization OPF problem. In the most recent literature 34 , Li and Gong formulate the optimal power flow with stochastic wind and solar energy as a multi-objective optimization problem and a multi-objective evolutionary algorithm based on non-dominated sorting with constraint handling technique is presented to solve it. In 4 , tidal energy is also added as a new renewable energy source in the optimization of multi objective OPF problems. In 35,36 , novel heuristic algorithms Slime Mould Algorithm and Hunger Games Algorithm were used to solve single and multi-objective optimal power flow problems. The summary of related studies is listed in Table 1. Although this literature has studied OPF problems, including renewable energy, it is only the initial stage of research on OPF including renewable energy. This is because most studies in the literature only study the OPF problem with a single optimization objective in the grid with renewable energy supply, or use the linear weighting method to integrate multiple objectives into one optimization objective. It's important to reduce emissions and other indicators while reducing costs to the economy. Therefore, it can be concluded that the multi-objective OPF problem with renewable energy is worthy of further study 37,38 .

Mathematical models
The OPF problem can be regarded as an optimization problem, classical OPF problems are single-objective. In contrast, the multi-objective OPF problem has many advantages in its solution. It can achieve the optimization of multiple optimization goals in one solution, such as fuel cost emission, power loss and voltage deviation. Multi-objective is not only a solution, it is a solution set, which can give decision-makers more opportunities to choose a compromise solution. The mathematical model of the multi-objective OPF problem can be defined by the following: where N is the number of objective functions, f i is the objective function to be optimized in a multi-objective OPF problem, i = 1,2,…,N; a and b are the control and state variable vectors, respectively. Multiple constraints need to be satisfied in the MOOPF problem, g (a, b) and h (a, b) represent the equality constraint and inequality constraint in the multi-objective OPF problem. The goal of solving the MOOPF problem is to find an optimal control variable a, which minimizes fuel cost, emission, power loss and voltage deviation. The Control variable vector is given in (3).
where P TG is the active power of the thermal generators; V TG is the voltage value of all generator unit buses; Q SH is described as the shunt VAR compensation. NG is identified as the number of thermal generator buses in the test network, and NC is identified as the number of shunt compensators in the test network.
The state variable vector is given in (4) (1) www.nature.com/scientificreports/ where P TG 1 represents the swing slack generating unit, V m r indicates the voltage magnitude at the r-th load bus, NL is the number's value of load buses; Q TG i is regarded as the reactive power outputs at the i-th generator bus, D L is the apparent power of the transmission lines, nl is the number's value of the transmission lines.
Constraints. The power system in the MOOPF problem needs to satisfy many power flow constraints, including equality constraints and inequality constraints.
Equality constraints. The equality constraint stems primarily from the reality that the generator's active power must equal the active load demand and active power loss. The equality constraints of MOOPF problem can be defined as: www.nature.com/scientificreports/ where P D i and Q D i represent the active powers and the reactive load demands of the i-th load buses. θ i define the i-th bus voltage angle. G ij and H ij are identified as the conductance and susceptance values of the transmission line between the i-th and j-th buses. ND is the number of buses, and i belong to bus number 1 to ND.
Inequality constraints. The inequality constraints of the MOOPF problem are described as: (a) Generator constraints: (b) Shunt compensator constraints: (c) Transformer constraints: (d) Contingency constraints: where T K regard as the k-th branch transformer tap. P TG i , Q TG i , V TG i , Q SH j , T K , V m r and D L n must limit between its upper ( P max The cost of energy spent in the system. Part of thermal generator set. Considering the valve point loading effect, the total cost of thermal power units ( C TG ) is calculated as follows: where C TG regard as the thermal generation cost, a i , b i , c i are the cost calculation coefficients generated by the i-th thermal generator set, d i , e i are the cost calculation coefficients of the i-th thermal generator set due to the valve point loading effect. The specific values of these parameters will be given in the experiment in 34 .
Part of direct cost of wind, PV and tidal. Wind, photovoltaic and tidal power generation require the purchase and installation of equipment, so operators must pay this cost, which becomes the direct cost of renewable energy. These costs are proportional to the power of the equipment and can be calculated from (8), (9) 29 and (10) 39 .
where α i , β j and l represents the purchase and installation cost coefficient of the i-th wind power turbine, the j-th Solar photovoltaic panel and the l-th tidal power plant. P ws,i , P pvs,j and P ts,l regard as the scheduled power of the i-th wind power farm, the j-th photovoltaic power plant and the l-th tidal power plant. www.nature.com/scientificreports/ Part of uncertainty cost of wind, PV and tidal. Due to the uncertainty and intermittency of wind power, photovoltaic power generation and tidal energy, there will be two situations in which the demand for electricity is greater than the supply of renewable energy and the demand for electricity is less than the supply of renewable energy. When the first situation occurs, the power dispatching system needs to raise other thermal generating units to supplement the power supply to reach a state of equilibrium, so this situation will generate additional supplementary generation costs. In the second case, there will be a surplus of electricity generated from renewable sources, so the power dispatch system will have to pay for the extra renewable power 29 . The expenses paid by operators in these two instances are reserve expenses and penalty expenses, respectively. Equation (11) can be used to determine the uncertainty cost of wind power 29 .
where C Wc is the total cost of wind energy source unknown. C Rw and C pw denotes the wind turbine's reserve and penalty costs, respectively. K Rw,i and K Pw,i indicate the i-th wind turbine's reserve and penalty cost coefficient, respectively. Nw denotes the total quantity of wind turbines in the power system. P wr,i is the rated output capacity of the i-th wind farm plant. P ws,i is defined as the planned out power from i-th wind farm plant. P w,i is the output capacity of the i-th wind farm plant. Uncertainty cost of PV power can be calculated by (12) 29 . where C Sc is the total cost of solar electricity supply uncertainty. C Rs and C Ps represents the reserve and penalty expenses of the solar energy, respectively. K Rs and K Ps represent the reserve and penalty cost coefficient of solar energy source, respectively. P s,i is represent the actual generating capacity of the i-th solar power plant. f s (P s,i > P pvs,i ) and f s (P s,i < P pvs,i ) are defined as the probability of actual generating capacity more and less than the scheduled power, respectively. E(P s,i > P pvs,i ) and E(P s,i < P pvs,i ) represent the PV power expectancy above and below the scheduled power of the i-th PV power plant, correspondingly. Uncertainty cost of tide power can be calculated by(13) 4,39 .
where C Tc is the total cost of tidal electricity supply uncertainty.C Rt and C Pt represents the reserve and penalty expenses of the tidal energy. K Rt and K Pt represent the reserve and penalty cost coefficient of tidal energy source. P t,i represent the actual generating capacity of the i-th tidal power plant. f t (P t > P ts ) and f t (P t < P ts ) are defined as the likelihood of actual generating capacity being greater than or less than the planned tidal power plant capacity, respectively. (P t,i > P ts ) and E(P t,i < P ts ) are the expectancy of tidal power plant above and below the scheduled power of i-th tidal power plant, respectively.
Objective function. The optimization of the MOOPF problem primarily includes economic and environmental optimization, and this paper employs four objective functions to measure and optimize economic and environmental indicators. The four objective functions are total generator cost, emission, power loss and voltage deviation. The total generators cost is defined as follows: In order to promote the use of green energy, regulators need to assess the pollution generated by thermal generators. Therefore, the second target to be optimized is the emissions of thermal units in the power grid, which can be evaluated by (15).
where, m i , n i , w i , t i and r i represents some emission coefficients in the i-th thermal generator.
Reducing the active power loss in the line is also an important optimization objective. The third objective to be optimized is the active power loss, which can be calculated as follows: www.nature.com/scientificreports/ Voltage deviation has a great influence on the voltage quality in the power system, so voltage deviation will be the fourth target to be optimized. Voltage deviation in the MOOPF problem can be calculated as follows: According to the description of the model, to solve the MOOPF problem, we need to search for a control variable vector a, and the make the four objective functions minimized at the same time. In addition, a needs to satisfy the constraints. Because solving the MOOPF problem with renewable energy is extremely difficult, we suggest a general framework for solving it. The general framework for solving the MOOPF problem using metaheuristics is shown in Fig. 1.

Uncertainty and power models for renewable energy
The stochastic modeling of renewable energy in the model will be described in this section.
Wind speed, solar radiation and tidal flow probability. Wind speed is characterized by intermittent and uncertainty. It is well known that probability density functions can be used to calculate the mean output of wind turbines [28][29][30] . The wind speed (v) m/s follows a Weibull probability distribution function and is calculated using the scale parameter ( γ ) and form parameter ( φ ) as shown below. Many researchers have studied the probability distribution of wind speed, and the Weibull distribution 28 is considered to be a probability model that can well fit the distribution of wind speed. Wind speed likelihood can be determined using the Weibull probability density function as follows: where v is defined as the wind speed, φ and γ are represents the shape and scale parameters. With the values provided in Table 5, φ and γ were carefully selected to ensure both diversity and realistic geographic locations for wind farm sites.
Create initial population randomly: www.nature.com/scientificreports/ The study discovered that where meteorological conditions are more dispersed, the lognormal function accurately represents the frequency distribution 40 . Use the lognormal probability density function to model solar irradiance 40 for the solar probability, which can be described as follows: where S represents the solar irradiance, δ and µ are equal to 0.6 and 6, respectively.
In earlier works, Gumbel distribution was typically used to calculate the probability model of flow rate in the tidal range 4,39 . Equation (20) is calculates the likelihood of tidal energy transfer.
From the result reported in 33 , the IEEE 30-bus system was modified to include the supply of wind energy, solar energy and tidal energy. In order to show the comparability of the results, the parameters of the probability distribution function from the previous literature are used in this paper. Different parameters can enhance the diversity and uncertainty of the renewable energy supply. The PDF parameters of these renewable energy sources are detailed in section "Result on the modified IEEE 30-bus with renewable energy".
Wind turbine and solar panel power model. Wind turbines and solar panels will provide electricity for the electrical network. Therefore, their power output needs to be calculated according to the wind speed probability and solar radiation probability in the previous section. The turbine actual output power is a function of wind speed 28 , the function can be formulated as follows: where v in = 3 m/s,v out = 25 m/s are defined as the cut-in and cut-out wind speeds, v r = 16 m/s represents the rated wind speed, P wr = 3 MW is the wind turbine rated output power. The turbine has three states in wind speed. In the first situation, v < v in or v > v out , the wind turbine will be stationary or locked to protect the speed does not exceed the limit rotor speed. In the second situation, v in ≤ v ≤ v r , the wind turbine will output power, according to wind speed. In the finally situation, v r < v ≤ v out , the wind turbine will be continuing to produce electricity at rated power. According to, the probabilities of three different cases can be calculated as follows: According to, the solar panel's electricity output is a function of solar irradiance (S), which is defined as follows 40 : where, P sr is the rated output power of the solar panel unit. S std = 800 W/m 2 is the standard environment's solar irradiance.R c = 120 W/m 2 is a certain irradiance point.
Tidal power, which generates electricity, was used by the sea water enters the reservoir at high tide to turn turbines, and then retreats from the reservoir back to the sea at low tide to turn turbines. Figure 2 shows the process of using tidal energy to generate electricity during high and low tides. The generating capacity of turbines in a tidal power plant can be calculated by (26) 39 .
where ρ is the water density (kg/m 3 ), g is the gravity acceleration (m/s 2 ), Q is the discharge value (m 3 /s) across the turbine set. ε is the turbine efficiency, H is the difference in height between the reservoir and the sea surface. These parameters of the tidal system sere set as H = 3.2 m, ρ = 1025 kg/m 3 , ε = 0.85 and g = 9.81 m/s 2 .

Multi objective pathfinder optimization algorithm (MOPFA)
The MOPFA will be suggested first in this part. Single Pathfinder Algorithm was proposed by Yapici and Cetinkaya in 2019, PFA has a strong global search capability, but the algorithm only deals with the single objective optimization problem. So in this research, we proposed the MOPFA, and then applied it to solve the MOOPF problem. Similar to the multi-objective particle swarm optimization algorithm (MOPSO) 41 . The first component is the archive, which is used to store the Pareto optimal solutions found so far. The pathfinder role selection process, which chooses the most likely pathfinder from the database and can guide the group to the optimum area, is the second component.

Pareto dominance.
When faced with solving a multi objective optimization problem 42 , x 1 and x 2 are two solutions, f i is the i-th objective function in this problem. If all the values of the objective function calculated by x 1 have the same or at least one better value than all the values of the objective function calculated by x 2 , x 1 is said to dominate x 2. This relationship can be symbolically expressed as x ≻y . If there is no x in feasible solution set makes x≻x ′ , x ′ is defined as Pareto optimal solution. All the Pareto optimal solutions are combined into a set, which is Pareto optimal solution set. In addition, the set of values of the multiple objective function corresponding to the Pareto optimal solution set is called the Pareto front 42 .
PFA population initialization. The MOPFA's first stage is to initialize the population so that it is evenly dispersed in the search area, which is accomplished using (27): where, x is the position vector of the individual population, LB and UB are the upper and lower bounds of the problem.
External archive initialization and update rules. The external repository holds the non-dominated optimal solutions as well as the collection of non-dominated optimal solutions found prior to the current run. The archive's capacity is fixed, typically half the size of the populace. The non-dominated optimal solutions derived from the initialize population were appended to the external archive component when it was started. The archive revised criteria can be specified as follows in subsequent iterations: In situation 4, crowding distance is a measure of the distance between one non-dominant solution and other adjacent non-dominant solutions in the archive. Equation (28) can be used to compute the crowding distance.
where, Max and Min represent the maximum and minimum values of each objective, respectively. Archive_size is the archive's capacity. According to (28), the crowing distance index d indicator is defined as the number of neighboring solutions that are less than distance d. A roulette technique was used to arbitrarily delete a solution from the external archive in order to add a new non-dominated solution to the complete external archive. Assign a (27) x = LB + rand(0, 1) · (UB − LB) www.nature.com/scientificreports/ probability to each non-dominated solution in the archive according to the crowding distance indicator, this fitness calculation is completed by (29).
where, P i is the probability of a non-dominant solution in an external archive is selected to be deleted by roulette method.N sum is the sum of the crowding distance index of each non-dominant solution in the external archive.
Pathfinder individual update rules. A non-dominant solution set stored in an external archive can be regarded as an elite individual, in the pathfinder algorithm, the population was divided into followers and pathfinders, and pathfinders led the population to the most promising region. Because elite archiving and the pathfinder's leader behavior are both elite leader behaviors, treat the external archive as the pathfinder individual in the PFA and update it using the pathfinder's individual update rules in the PFA and the optimal protection strategy. The pathfinder updating position is obtained from (30).
where x i+1 p indicates the place of the i + 1 generation pathfinder, x i p is the location of the i-th pathfinder. x i−1 p the i − 1 generation pathfinder's position, i represents the number of current iterations, r 3 represents a random integer from a uniform distribution at [0, 1]. η is derived from (21).
where i max indicates the highest number of iterations,u 2 is a random integer evenly spread between [− 1, 1].
Follower individual update rules. The external archive stores the non-dominant solutions found so far, which can be regarded as people in the population's most hopeful area. The pathfinder individual description in Single objective PFA is matches. Different from PFA, in MOPFA treat population individuals as followers and external archives as pathfinders. According to the size of follower population, the corresponding pathfinder matrix is constructed, follower's update process as follows: where i represent the current iteration, x n represents the position of the follower of the population. W 1 , W 2 are two randomly generated vectors calculated used Eqs. (33) and (34), r 1 and r 2 are a uniformly distributed random number generated randomly between [0, 1]. W 1 ,W 2 can control the weight of the follower moving to the pathfinder and the neighboring individuals in the population. ε is the vibrancy vector, and its calculation can be obtained from (35), ij is the distance between the i-th and the j-th position in population.
The pseudo code of MOPFA is given in Algorithm 1, and the flowchart of the multi-objective pathfinder algorithm is shown in Fig. 3.

Experimental results and analysis
In this part, MOPFA is used to answer the multi objective optimal power flow issue (MOOPF), and the experimental findings are thoroughly examined. Generation cost, power loss, voltage deviation and emission were studied as objective functions in this research. In order to study the emission reduction effect of renewable energy in the MOOPF problem, renewable energy is included in the energy source of the power network. The experiment was conducted on an adapted IEEE 30-bus. The power flow formulae of the suggested MOOPF model with renewable source were calculated using MATPOWER 6.0. Four cases were studied for different optimization purposes. Case 1 is to reduce both production costs and emissions at the same time. Case 2 and Case 3 are focused on three objectives. The most difficult is Case 4, which aims to optimize all four objective functions at the same time. The specific experiments case content is listed in Table 2. The optimization results of MOPFA and several other multi-objective optimization algorithms on modified IEEE 30-bus for MOOPF problem with four cases are listed in this section. The Hyper-Volume (HV) indicator 43 obtained by MOPFA and other multi-objective optimization algorithms under four optimization objective cases were compared and also the results were also compared using statistical test analysis. Using the well-known fuzzy decision system 44

HV indicators.
In the multi-objective optimization algorithm, the quality of the algorithm needs to be measured through a variety of indicators. The HV indicator is a comprehensive evaluation indicator, and can be compared without the real Pareto frontier. The higher the HV value is, the better the algorithm effect will be.

Fuzzy decision system.
In this paper, the following fuzzy decision system is used to select the compromise solution in the Pareto optimal solution set. The calculation method is as follows 44 : where γ k n is the membership function value of n-th objective for k-th non-dominated solution; f k n is the fitness value of n-th objective for the k-th non-dominated solution; f min

Result on the modified IEEE 30-bus with renewable energy. To integrate renewable energy into
the grid's electricity supply, IEEE 30-bus was modified in this study. According to Ref. 28 , the system consisted of 41 transmission lines, 6 generating units, 9 shunt VAR compensators, and 4 transformer tap settings and its total active and reactive load demands were 283.4 MW and 126.2 MVAR, respectively. In the modified IEEE 30-bus the thermal generators at buses 5, 11 were replaced by wind generators respectively, wind generator farm www.nature.com/scientificreports/  www.nature.com/scientificreports/ at bus 5 has 15 turbines, and bus 11 has 10 turbines. The thermal generators at buses 13 and 8 were replaced by solar generator and tidal generators respectively, tidal generator at bus 8 has 4 generating sets. Table 3 describes the modified IEEE 30-bus settings, and the renewable uncertainty cost coefficients are listed in Table 4 29 . The modified IEEE 30-bus structure diagram was shown in Fig. 4. Parameter values of the probability distribution function of simulated renewable energy in all Cases 1-4 are given in Table 5.
Result on case 1: minimize the generator cost and emission. In this case, the goal of optimization is to minimize generation costs and emissions. Cost and emission are in conflict, but power system generation cost and emission control should also be given more attention to. After the 30 runs, the Hyper-Volume (HV) indicator obtained by MOPFA and other algorithms in Case 1 are listed in Table 6, which statistically analyzes the maximum (Max), average (Mean) and minimum (Min) values of each algorithm in the 30 runs and the best results are highlighted in boldface. Furthermore, in order to better test the performance of the algorithm, the Wilkerson rank sum test is used to rank the algorithm. From Table 6, it can be calculated that MOPFA obtained the maximum HV indicator of 0.18025 and the minimum value of 0.1423 after 30 runs, which was better than other algorithms. The average value was also compared, MOPFA achieved 0.14773, ranking first. Although the mean square deviation value was not as good as MOAHA's, the gap was also small. MOPFA's Wilkerson rank sum test score is 4.87, ranking it first among the six multi-objective algorithms. The boxplot figure of the HV indicator obtained by each algorithm on Case 1 is given in Fig. 5. The Pareto front obtained by the HV index maximum runtime algorithm is plotted in Fig. 6. It  www.nature.com/scientificreports/ can be seen that the Pareto front obtained by MOPFA is more uniform and MOPFA can give a more diverse set of Pareto optimal solutions. According to the fuzzy decision system mentioned in the preceding part of this paper, the compromise solutions were selected from the Pareto front in Fig. 6. The compromise solutions on Case 1 of MOPFA and other well-known multi-objective optimization algorithms are reported in Table 7. The composition of renewable generation energy and thermal generation energy in the solution obtained by each algorithm is shown in Fig. 7. The compromise solution obtained by MOPFA is marked in Fig. 6, the value of generation cost is 813.8379 ($/h) and emission is 0.2057 (t/h). From Table 7, it can be calculated that MOPFA obtained the best emission value: 0.2057, but MOMVO obtained the best generation cost value: 800.0093. From Fig. 7, it can be seen why MOPFA's solution did  www.nature.com/scientificreports/ not achieve the best cost. Because MOPFA's solution uses the total thermal power and renewable energy is 139.1494 MW ( P G 1 + P G 2 ) and 148.3646 MW ( P G 8 to P G 13 ), while MOMVO's solution uses the total thermal power and renewable energy is 141.5587 MW ( P G 1 + P G 2 ) and 143.4163 MW ( P G 8 to P G 13 ). MOPFA's solution uses more renewable energy and renewable energy uncertainty will increase the generation cost,  Table 3. Compare MOPFA's compromise solution with the published literature was reported in Table 8, MOPFA's generation cost is lower than other algorithms except ACNSDE-SF 30 , but MOPFA's emission is lower than ACNSDE-SF. Both generation cost and emission are considered, MOPFA obtained value is lower than MODFA 32 , MOEA/D-SF 33 , ESDE 19 , PSO-SSO 34 and MOMICA 35 . In summary, the solutions given by each algorithm are not dominated by each other, and the decision maker decides which one to choose, but MOPFA's HV index of the Pareto front is higher, indicating that the Pareto front is more evenly distributed, and it can provide a greater diversity of solutions. Moreover, MOPFA's compromise solution is more inclined toward the utilization of renewable energy, which will reduce the spending of enterprises when the government imposes higher emission taxes.
Result on Case 2 and Case 3. In Table 9, for Cases 2 and 3, the number of optimization objectives is three. Case 2's goal is to reduce generation costs, emissions, and real power loss. Case 3 aims to minimize generation cost, emission and voltage deviation. In Case 2 and Case 3, MOPFA obtained the maximum values of HV index 0.108043 and 0.133712 in the six algorithms and also achieved the best results on the minimum and average values of the HV indicators. The Wilkerson rank sum test scores of MOPFA in Cases 2 and 3 were 5.40 and 5.93, respectively, MOPFA ranked first among six algorithms in the two cases. The boxplot figure of the HV indicator obtained by each algorithm in cases 2 and 3 is given in Figs. 8 and 9, the variance of MOPFA is very small, indicating that the solution obtained by MOPFA is very stable. The Pareto front, resulting from the maximum run of the HV indicator is plotted in Figs. 10 and 11, it is obvious that the Pareto front obtained by MOPFA in Cases 2 and 3 is more evenly distributed than that obtained by other algorithms. The compromise solutions obtained by each algorithm in the two cases are given in Tables 10 and 12, respectively. In Case 2, the compromise solution obtained by MOPFA is marked in Fig. 10, the value of generation   www.nature.com/scientificreports/ cost is 863.2328 ($/h), emission is 0.0508 (t/h) and real power loss is 2.3446. In Case 3, the compromise solution obtained by MOPFA is marked in Fig. 11, the value of the generation cost is 884.6613 ($/h), emission is 0.0489 (t/h) and the voltage deviation is 0.4037. Figures 12 and 13 provided the composition of renewable generation energy and thermal generation energy in the solutions obtained by each algorithm. From Table 10, it can be seen that MOPFA achieved the best real power loss of 2.3446 MW, while the second-best score of 2.9515 MW was obtained by MOAHA, MOPFA's solution results in a 30% reduction in power loss compared to MOAHA's solution. MOSSA obtained the best generation cost value 834.2038 ($/h) in Case 2. MOAHA obtained the best emission value 0.0482 (ton/h) in Case 2. It can be seen from Fig. 11, MOPFA uses 81.9 MW thermal power generation energy, which is much lower than the 127.4 MW, 111.2 MW, 125.9 MW, 105.3 MW and 138.2 MW of other algorithms. In particular, wind farms, photovoltaic plants and tidal power plants connected to buses 5, 8, 11 and 13 are scheduled to be assigned 74.5949 MW, 57.8024 MW and 47.2424 MW respectively. This is more than 90% of the intended capacity of these renewable power plants, which will undoubtedly increase the default and reserve costs, but as the experimental data shows, the renewable energy generation process does not require fuel costs, so the cost increase is acceptable. The Case 2 results obtained by MOPFA were compared with the recently published literature in Table 11. The comparison results show that MOPFA obtained the minimum emission value 0.0508 (ton/h) than other algorithms. Although MOPFA's power loss value of 2.1891 (ton/h) is inferior to TLBO's 38 , the cost and emission of MOPFA are significantly lower than TLBO's 38 among the three targets. MOPFA's emission and power loss is lower than NASGA-II-SF 30 , MOEA/D-SF 33 , PSO-SSO 34 , MOAGDE 31 and ACNSDE-SF 30 . Different from Cases 2 and 3 adopts voltage deviation instead of power loss optimization target. From Table 12, it can be seen that MOPFA provides the Table 9. HV-indicator in Case 2 and Case 3. Significant values are in bold.   Table 13. MOPFA achieves the best maximum, minimum and average values, the HV indicator variance plots of each algorithm are plotted in Fig. 14. MOPFA's results are remarkably stable and better than those of the other five algorithms, the results of the 30 runs were ranked by Wilkerson rank sum and MOPFA ranked first in six algorithms, it's score is 5.88. Figure 15 shows the best Pareto front obtained by each algorithm. From Fig. 15, it can be seen that the distribution of Pareto optimal solutions of MOPFA on each objective is very uniform,   Table 14. According to this table, MOPFA's solution has the lowest emission: 0.0486 (ton/h) and the highest real power loss (3.0052). MOMVO obtained the minimum value for generation cost, but the value of emission and real power loss are more than twice as high as MOPFA's. Through the energy composition histogram in Fig. 16, the reason for the compromise solution is analyzed. Because MOPFA's solution dispatches the most tidal energy, this may be because tidal power is more stable than wind power, reducing the extra cost of renewables to some extent. The total renewable power of MOPFA's solution is 178.9 MW, it is bigger than NASGA-II: 147.9 MW. But dispatching    www.nature.com/scientificreports/ through MOPFA increases the use of more stable tidal energy in solutions, which reduces the uncertainty of wind and photovoltaic energy. As a result, MOPFA's generation costs 5.6% less than NASGA-II. Table 15 reports a comparison of the solutions obtained by MOPFA and the recently published literature. From Table 15, it is obvious that MOPFA gets the smallest emission (0.0486), PSO-SSO 50 gets the smallest cost, NASGA-II-SF 34 gets the best voltage deviation. The solutions provided by these algorithms do not dominate each other, but the MOPFA solution should be considered when considering emissions taxation. In summary, MOPFA offers a more diverse set of solutions, and MOFPA offers solutions to increase the use of more stable tidal energy when renewable energy uncertainty increases the cost of the solution. This will help to resolve the contradiction between costs and emissions, and further demonstrates that the use of three renewable energy sources connected to the grid can greatly reduce emissions, and this energy composition will be very promising in the future.
System constraints of the best compromise solution. In the MOOPF problem, the voltage on each branch must be within the secure range to guarantee the power system's regular functioning. Because the bus containing the green energy units and thermal power units is already in the range, the voltage on the remaining 24 buses must be in the range of [0.95 p.u., 1.05 p.u.]. As a result, in order to validate the effectiveness of the answer achieved by each algorithm in Cases 1-4, Fig. 17 depicts the voltage of each branch of each algorithm on the IEEE 30-bus in the four study cases, clearly showing:  www.nature.com/scientificreports/ In Cases 1-4, the intermediate solution produced by six algorithms meets the branch voltage restriction, demonstrating the possibility of using a meta-heuristic algorithm to solve the MOOPF issue. When we concentrate on Case 1, each algorithm's obtained voltage in bus 3 is near to, but not exceeding, the upper bound. The voltage intervals of each algorithm are not noticeably different in the other three testing instances. In summary, each calculation's compromise answer meets the branch voltage limit.

Conclusion and future work
In this work, a multi-objective optimization power flow (MOOPF) problem with stochastic wind, solar power and tidal power models is introduced. Different scenarios of renewable energy supply are simulated by considering different probability distribution functions. To solve this complex multi-objective optimization problem, a novel multi-objective optimization algorithm MOPFA was proposed. The thermal power generators in IEEE 30-bus were replaced by wind turbines, photovoltaic power plants, and tidal power generation equipment. The simulation is carried out on the modified IEEE 30-bus system. The experiment was conducted in four different cases and the experimental results are compared with those of other well-known multi-objective optimization algorithms. Statistical results showed that MOPFA achieved the best HV indicator in all four cases. In addition, MOPFA's Wilkerson rank sum test was also ranked first, while MOPFA is slightly better than other algorithms. In solve the multi-objective optimal power flow problem, MOPFA can obtain a more widely distributed solution set, and the solution that satisfies the constraints, so it can be used as the preferred algorithm to solve this problem. The compromise solution is calculated from the solution set obtained by MOPFA by fuzzy logic, and the compromise solution of MOPFA uses more renewable energy supply and effectively reduces emissions. Incorporating renewable energy into the power system can reduce emissions while maintaining system stability,  www.nature.com/scientificreports/ which will be an advantage in the future, so the compromise solution obtained by MOPFA is more in line with the development trend of future power dispatch. In the future, this research will focus on the direction of solving the MOOPF problem on larger IEEE test systems with real wind energy datasets, and search for real power price datasets from the government to build more accurate models. The composition of renewable and thermal power in case 4 www.nature.com/scientificreports/  www.nature.com/scientificreports/

Data availability
The datasets used and analysed during the current study available from the corresponding author on reasonable request.